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MODEL 



JEAN BOURGAIN* 

Abstract. The purpose of this Note is to provide a deterministic implementa- 
tion of the random wave model for the number of nodal domains in the context 
of the two-dimensional torus. The approach is based on recent work due to 
Nazarov and Sodin and arithmetical properties of lattice points on circles. 



1. Introduction 

This Note originates from the work of Nazarov and Sodin ( |N-S| and [S]) on 
the behavior of the number of nodal domains of random eigenfunctions at high 
energy. It was sown in |N-S] that the number Ne of a random eigenfunction os 
of eigenvalue E obeys the so-called random wave model (RWM) for large E and, 
with large probability, the ratio 47r^ is close to a constant cr > 0. According to 
the Bogomolny-Schmit jB-Sl| . |B-S2| prediction, this number a can be computed 
based on a bond percolation model leading to a conjectured value 

a ^ ^^'^ ^ 0.0624... (1.1) 

TT 

While the IN-S| work establishes in particular the positivity of a, it does not shed 

any light on its actual value. Note that (II. ip is considerably smaller than the 

general (deterministic) upper bound provided by Pleijel's inequality 

N /2\2 ^ ^ 

limsup— <(-] =0.691.. (1.2) 

n— >C30 1^ J 

with j the smallest positive zero of the Bessel function Jq and n x the wave 
number (see also ^ for a small improvement). 

Better upper bounds on a may be obtained by evaluation of certain geometric 
parameters using Kac-Rice type arguments. It was shown in particular in [K] that 
(T < -i- = 0.225... by computation of the expected number of horizontal tangencies 
to the nodal set. The same bound may be gotten from its expected total curvature 
(cf. [Bi^ l. 

In what follows, we do not intend to study further the RWM or the Bogomolny- 
Schmit heuristics. Rather, we are interested in a deterministic implementation of 
the RWM in certain situations. The idea is very simple. Assume —A/ = Ef 
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an eigenfunction for large E. Fixing some base point x G M in the manifold E, 
we are considering restrictions fx of / to neighborhoods of x of the order 0(-^) 
(in fact with R slowing growing to infinity with E). In certain instances, one 
may then show that the ensemble {fx)xeM resembles that of a Gaussian random 
wave function. It turns out that for M — the 2-dimensional flat torus and 
eigenfunctions 

fix)^ J2 flce(x.e) (e(a)=e2"-) (1.3) 

(a_^ — a^) where £e — {£, Z^; — E} satisfies suitable arithmetical assump- 
tions and X]|^|2=_E I^cP'^c/x/b becomes well-distributed on the unit circle, this idea 
may be worked out rather easily. On the arithmetic side, we rely on angular equidis- 
tribution results |E-H] (see also jF-K-W| and related references) and also the recent 
work (B-B| on additive relations in the set £e- Naturally, one runs into stability 
problems for the number of nodal domains when perturbing slightly the eigenfunc- 
tions, but these analytical issues have been already completely addressed as part 
of the remarkable work of Sodin and Nazarov. In particular, extensive use is made 
from the results in [S]. 

Recall also that from a result due to A. Stern [St] (see also [D), there is no 
nontrivial lower bound on the number of nodal domains for i? — > oo, which may 
equal two. Thus for eigenfunctions (|1.3p . some further assumptions are needed. 
Possibly, the equidistribution of the measures J2\^\2=e \^i\'^^^/\/E '-'^ suffice, 
but we are only able to establish this in certain cases (for instance assuming E has 
a bounded number of prime factors and also in a statistical sense, i.e. for 'most' 
E). 

Beyond the arithmetical input and the results from [S] , our analysis is essentially 
straightforward. No effort has been made to obtain quantitatively more refined 
results. A more general outlook on the approach is discussed in the last section. 

Let us return to our model and be more specific. 

Assume E G Z+ a large odd integer which is a sum of 2 squares; we assume 
moreover E of the form 

E^YIpI- (ea>l) (1.4) 
where its prime factors pa = l(mod4). Denote 

S = £e = {^ = (^1,6) e Z"; + ^2 = E}. (1.5) 

Identifying (^i, ^2) with the Gaussian integer £,i + i£,2 (z Z -\- iL and denoting pa — 
TTaTTa the factorization of Pa in Gaussian primes, the set £ is obtained as 
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up to multiplication by ±1 and ±z. In particular 

\£\^4.Y[il + e^) = W. (1.7) 



Writing E = \^ and 



we obtain 



whit angles 



^a = ka|e'^° (1.8) 

C = Ae*'^ for ^^£e (1.9) 

Tp = ^(2ja - ea)9a and < ja < Cq, (1.10) 



up to multiples of f . 

The eigenfunctions of the Laplacian — A on with eigenvalue E are obtained 
as trigonometric polynomials 

/=^a4e(x.O. (1.11) 
Let us consider for simplicity the eigenfunction 

Our considerations in the remainder of the paper carry over verbatim to the sit- 
uation (|l.lip with |a^|,^ G £ equal and more general statements with arbitrary 
coefficients (aj)jg£ will be discussed later. 

Our aim is to show that under suitable assumptions on i? — > cx), the number of 
nodal domains of (I1.12p obeys the RWM. These assumptions are of arithmetical 
nature and may be loosely formulated as follows 

(D) The points {A^^^;^ G £e} become equidistributed on the unit circle for 
E ^ oo. 

(I) There are not to many non-trivial additive relations among the elements of 
£. 

While we only need (D) without further quantification, a more precise form of 
(I) will be required (See Definition [T] and Proposition [1} . Properties (D) and (I) 
may be addressed by classical results in number theory. By (|1.10p . (D) relates to 
angular distribution of Gaussian primes and we will refer to the results from |E-H) . 
A powerful tool to deduce bounds on the number of additive relations is provided 
by [E-S-S| on unit equations 

aia + ---+a«C« = 1 (1-13) 

with ^1, . . . , taken in a multiplicative subgroup G of C* of bounded rank (though 
the available results require some further assumption on the number of prime factors 
of E to be applicable to our problem). Alternatively, one may use the 'statistical' 
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results on additive relations proven in |B-B| to treat the case of 'typical' E. Precise 
statements will be given in section 4 (Theorems [21 [3l |4]). 

2. Local Analysis of the eigenfunction 
Let be equipped with normalized measure and let 



/(x) = ^^e(x.O (2.1) 



with £ = £e, E ~ and W — \£\. We always assume W oo with — > oo. 

In what follows, we will need several parameters, chosen in a particular order, 
that will be viewed as 0(1) for fixed E and eventually will tend to infinity with 
_B — > oo at sufficiently slow rate. 

Let \ <^ K = o{W) be a first large parameter and subdivide AS*^ in arcs of size 
leading to a corresponding partition 

K 

E^\jE'^^\ (2.2) 

fc=i 

More specifically, we subdivide the first quadrant of AS*^ and partition the other 
regions by reflection and symmetry. According to (D), assume that 

l--e,)w <\E^^^<[^+e,)w (2.3) 

for each k — \,...,K, where £i = ei{K). Choose a point ^^'^^ € letting 
^(fe') ^ _^(k) if £(/c') ^ 

Let i? 3> 1 be another parameter and denote 

(■{k)^R^ik)^ (2.4) 
A 

Hence |C^''''| = R- Fixing x £ T^, translate x by with y = (2/1,2/2) G [-5, \f 
and write 

i^.(2/) = /(x+f2/) =^E/fe(x,2/)e(C('=).y) (2.5) 

with 



Denote further 



^(2/) = = (2-7) 

fc=l 

with 

Cfe(x) =A.(x,0). (2.8) 
Our next goal is to show the following 
(i) For most x, Lpx is a perturbation of Fx considered as function of y. 
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(ii) The random vector {ck{x)}i<k<K with x ranging in has approximatively 
the same distribution as the Gaussian vector {gk}i<k<K , with gi, . . . ,gk 
IID normahzed complex Gaussians, subject to the reahty condition gk' = gk 
for ' = — (■('"^ At this point, we wih then be able to rely on the results 
from fS). 

Note that in (ii) , we should see K as fixed and the distributional approximation 
sufficient in order for the relevant Gaussian estimates from [S^ to carry over. 

Lemma 1. For any fixed s> 1, 

\\F^-vAlic^,<R'''K-\ (2.9) 

Proof. Since from standard Sobolev estimates, we may bound the C*-norm by the 
i/*+^-norm, it suffices to estimate 



K 

by I^M''^- □ 



It follows from (|2.9p that after fixing i?, we may ensure, taking K sufficiently 
large, that 

||i^.-^.||cs =o(l) (2.10) 

for most X G T^. 

We now turn our attention to the joint distribution of the vector {Ck{x)]l < 
k < K}, 

Ck{x) J2 ^(^■^) (2.11) 

when X ranges in T^. 

Switching notation a bit, it will be convenient to replace K by 2K and enumerate 

with Obviously c_fc = Cfc. 

We specify assumption (I) as follows. 

Definition 1. Fix < 7 < ^ and some B G Z+. We say that £ satisfies property 
1(^,3) if for 2 < £ < B , the number of non- degenerate additive relations of the 
form 

a + ---+e£=0 (2.12) 

among elements ^i, . . . , G £ is at most N^^ . By 'non- degenerate', we mean that 
in ()2.12p there is no proper vanishing sub-sum. 



There are various ways to select energies E for which above independence prop- 
erty holds and this will be addressed in a later section. 
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Definition 2. Let s > be a small parameter. Say that the random vector 
(ci,...,Cfe) where the Cj are <C- valued functions, \\cj\\2 ^ 1, is e-Gaussian, pro- 
vided for any (possibly unbounded) intervals /i, Ji, . . . , //f , Jk C M, we have 

mes[cieIixJi,...,CkelKxJK]--77rW / e~5(^?+y?+--+^?c+ai)d2;dy 

(2.13) 

Choosing e sufficiently small (in particular depending on K), (|2.13p will enable 
to approximate for 'nice' open sets fi C , mes[(ci, . . . ,ck) G ^] by the corre- 
sponding Gaussian measure. We prove 

Lemma 2. Given e > 0, there is B = B[K,e) such that if £ satisfies I{'j,B), then 
the vector function (ci, . . . , ck) on as defined above, is e-Gaussian. 

Proof. Well-known arguments reduce the problem to evaluating moments 

/ c^c^---c^c^ (2.14) 

with ri, Si, . . . ,rx , sk G U {0} and 

ri + si + • • • + TK + SK < Bi{K, e). 
Recall that the procedure consists indeed in evaluating the characteristic function 

/ g[Q!i Rgci+/3i ImciH hQ/i- Rcc^+^fc ^^'^'^^dx (2 15) 

with ai, Pi, . . . ,aKT Pk G IR, subject to some bound B2{k,e). Those arise by 
suitable truncations of the Fourier transform of intervals. Then, imposing some 
bound on |ci|, . . . , |cx|, Taylor expansion of the exponentials in (j2.15p leads to the 
expressions ()2.14p . 

Substituting ((2TT|) in ((2l4)) gives 

/ K \ h(ri+...+rK+si-\ \-sk) 

where (2.16) stands for the number of relations 

6,1+- • -Hi^r.-^'is S.'i,s,+- ■ -+^,1+- ■ •+avK-eki ^ks. - (2-17) 

with ^1,1,- •• ,Ci,ri and & ■ ■ ■ 

Trivial solutions to (|2.17p are those for which the multi-sets (i.e. taking into 
account multiplicities) 

{6,1) • ■ • )6,ri> ■ • ■ ,^K,l, ■ ■ ■ ,£,K,rK} (2-18) 

and 

coincide. Of course, to have a trivial solution requires 

ri = si, . . . ,rx = SK- (2-19) 



< e. 
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Otherwise, we call the solution non-trivial. 

Consider first the contribution of trivial solutions, assuming ()2.19p . 
Denote = T'^ and define ci, . . . ,ck on by 

= E e(V'«) with = (2.20) 

Recalling (|2.3p and taking into account the central limit theorem, the distribution of 
{ci, . . . ,ck) is approximatively Gaussian. The trivial solutions to (|2.17p contribute 
for 

' |5i|2'^^---|£Kr-= / \9ir'---\9kr''. (2.21) 



Consider next the contribution of non-trivial relations, which will be evaluated using 
our arithmetical assumption. Their number is obviously bounded by the number 
of nontrivial relations 

CiH 1-6 = with ^ = ri -)- \-rK + si-\ \- sk (2.22) 

in elements ^ from £. Partitioning (20) in minimal vanishing sub-sums, at least 
one of these relations will be non-trivial and therefore of length £' > 3. Property 
I{'-f,B), B > £, clearly implies the following bound 

C{e) J2 WW^^'^-'^"^ <C{(.)W^.W'^^^'"'\ (2.23) 

2iy<f-3 
i_ 

Multiplying with (■^) ^ , the resulting contribution of the non-trivial relations (|2.17p 
in (|2.14p is therefore at most 

B3(ii:,£).W^-(3-7) (2.24) 

which can be made arbitrarily small for W large enough. 

This proves Lemma [2] □ 



3. The number of nodal domains 
Consider the eigenfunction p.ip on T^. 

From general theory, the total length of the zero set Z{f) = {a; e T^; f{x) = 0} 
is 0{\) while each nodal domain has area at least 0(A^^). In particular, it follows 
that the number of nodal domains of diameter at least e^^A~^ is at most 0(e2A^). 
Here £2 is a small fixed constant. 

Choosing R sufficiently large, it clearly follows from the preceding that 




where N f is the number of nodal domains of / and Nf{x, p) the number of nodal 
domains contained in the open box x + (] — f,f[x] — f,f[). 
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Using our notation p.5|) . te first term on the rhs of p.ip equals 

^ / NfJx (3.2) 



with Np the number of components of Z{F) contained in ] — i, i[x] — |[- Note 
also that 

NF^=Nf[x,j) <OiR'). (3.3) 
Let (fi be defined by (|2.7p . According to Lemma [T] 

j \\F^-'pAc^ dx<R^K-\ (3.4) 
Hence, fixing another parameter £3 > 0, it follows that 

ll^x - ^xWci < £3 (3.5) 

except for j: in a subset F C of measure at most e^^R'^ K^^ < 82, taking K 
sufficiently large. Since then, by (13. 3p 

we may replace p.2p . up to 0(e2A2), by 

R^ Jy2 
where H-^xUci < ^3- 

For X G V, set ip^ — 0. Since the function ipx on satisfies — A(^a; = R'^(px, it 
follows again from the Faber-Krahn inequality that each nodal domain of ip,j. is of 
area at least 0{^) and hence N^^ < 0(i?2). Thus 



■ / NfJx < Oi\^.\V\) < esA^ 
Jv 



^Vx+^Ax dx (3.6) 



^ j^N^Jx<0{\^\V\)<e2\^ 



and in p.6p . the integral may be extended to T^. Consequently, we obtain 

A^ f 

~R^ y^^ (^■'^) 

The next step consists in invoking Lemma [21 which asserts that for W sufficiently 
large, the ensemble {^Px)xeT^ has approximately the same distribution as the Gauss- 
ian random function 

'^.^^Y.9u{io)e{e^.y) (3.8) 
with {gk\ IID normalized complex Gaussians subject to the condition 17^/ = for 

We claim that by choosing e small enough in Lemma [H one can replace p.7p by 

A2 



i?2 



iV<i>^+M,^da; (3.9) 
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where "^ui is some perturbing function, satisfying 

ll^^llci < 263 (3.10) 

and 

- A($^ + *„) = i?2($^ + $^). (3.11) 

Proof of the claim 

Take M ^ -^log K^l log — and subdivide the M-cube centered at in in 



cubes Qa of size £4 — Their number is 0[{AIe^ ^Y^)- 

For each a, denote 

Aa = {x e T^; (cfc(a;))^^^^^ e Qa} 

and 

with the probabihty space on which is defined. According to Lemma [U we 
can ensure that 

I I^qI — \Ba\ \ < £ for each a. 
Note that \Ba\ > 6{K, M, 84) and hence, for e smaU enough, we may ensure 

\\Aa\-\Ba\\<^\Ba\. (3.12) 

This permits to introduce subsets A'^ C Aa , B'^ C , such that 

\A'a\ = \B'J>{l-^)\Ba\ (3.13) 

and a measure preserving map 

Ta '■ B^ — > A^. 

Define on \JB'„ 

1 ^ 

^SfUy) = ^5][cfe(T„(c.)) -5,(c.)]e(C('=).2/) + V.„M (3.14) 
fe=i 

and set 

*„ = Oifw^UB;,. 

With this construction, 

duj 



= J2[ N^.+^Jx + 0{R'\{UB'J^\) 

„, A' 
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N^^+^Jx + 0{R\UA'J'=\) + 0{R'\(UB'^r\) (3.15) 
where we have used again that 

N^^+^^ <0{R^) and iV*„ <0{R^). 



Next 



\{^B'^Y\=Y^\B^\B'^\ 



+ |{^;^max^; l5fcMI > M}\ 

(3^3) £2 V- I R I , £2 £2 



and 



\{yjA'^r\ = + e T^; ^max^ |c,(a:)| > M}|. 

Again by ([333]) 

a 

From Lemma [5] 

if 

\{x e T^; max |cfe(x)| > M}\ < E mes[|cfe| > M] 

l<k<K 

<K{e+ mes[|gfc| >M]) < ^. 

Substituting in (I3.15P gives 

j N^^+^,Juj = j ^N^^+^Jx + 0{e2). (3.16) 
Finally, note that on , by p. 141) and choice of £4 

llf^llci < i?\/^£4 + IIV-t^mIIci < 2£3. 

Also, since either = or $^ + ^1';^ = (^2; + V'x for some 2;, it follows that 
-A($^ + 5'^) = R^i'^u + *tu)- This proves the claim. □ 

At this stage, we are reduced to study the expected number of nodal domains in 
] — ^,|[x] — i,i[of the perturbed Gaussian vector 

We make use of the work of Nazarov-Sodin and more specifically, several results 
from [S]. 

First there is the stability issue. Considering the random Gaussian function '^^ 
given by (|3.8p . clearly 

E^[|l<i>„|lc2] < 0(i?2). (3.17) 
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Invoking Lemma 5 from S , which is based on the independence of $aj and V$tj, 
we get some /? = 13{R, £2) > such that 

min max(|$„(x)|,|V$^(x)|) >/3 (3.18) 
for aU u! outside a set of measure less than hence contributing to 

for at most 0(62) ■ 

Property (13.18^ is crucial to derive a stability property for the number of nodal 
domains under perturbation (see [SJ, Lemma 6). Recall that the perturbation 
satisfies ||5'cj||ci < £3- Taking 



£3 =/3(i?,£2) 



1 



lOR 



(3.19) 



Lemma 7 from [S] applied with / = "Scjiff = and a — 26^ implies in particular 
the following 

Af$^+*^ > number of components of Z(^i^) contained in the square 
1 1 



Q 



2' 2 



2 1 2a 

and at distance at least — — > — from dQ. (3.20) 
2R p 



Note that since < £3, (|3.18l) also implies that 



/3 



minmax(|($^ + ^'„)(x)|,|V($<^+*c.)(a;)|) > ^ 



(3.21) 



Another apphcation of [5], Lemma 7 taking / = $0; + ^'^, .9 = — ^'w yields conversely 
that 



Number of components of Z(<^^) contained in 



1 1 1 



1 



2 2i?' 2 2R 
It follows from the two-sided inequalities (|3.20p . (I3.22p that 



E 



^ components C of Z[(j)^) contained in 



1 11 1 n2 

2^ 2R'2^ 2R. 



dist(C, dQ) < - 
K 

+ 0(£2). 



Recalhng dHSl) and Cfe = f -^fc with 6 e Z^, = A, we have 

R 



(3.22) 



for which 

(3.23) 



defining 



K 



fu>{x) = ^ Vgfc(a;)e(Cfc.x) 
fc=i 



(3.24) 
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and iV(0, y ; fuj) the number of components of Z{f^) contained in 

Qii ~] ~ 5I' 5x[^] ~ ll' lit- Thus fu, is a Gaussian random eigenfunction of of 

eigenvalue E. 

The first term in (|3.23p accounts for the number of components C of contained 
in Qr+1 and such that dist(C, SQi^) < j. Each of these components has area at 
least 0( jt) and it follows from the Kac-Rice formula that 

[length (Zf^ n Qr+i)] < o(^). (3.25) 



li 3 

# components C of Z{f^) s.t. C c dist(C, 5Q_r) < t- < 0{R^ logi?) 

A J 



From these facts, one easily derives that 
E 

. . 

(3.26) 

(we first exclude those components C of size at least logi?.^^, assuming 

length {Zf^ n Qb+i) < log-R.-^ ^^id then exploit the area lower bound of each 

component). 

Hence we proved that 

E[7V$^+*J =E[7v(0,p/^)] +0(e2 + i?'/'logi?). (3.27) 

The expectation of Af(0, y , /(^) in the limfl^.oo limA_^oo is given by Theorem 5 in 
[5] and we get in our situation 

R^.vip) (3.28) 

where v{p) is the constant given by [S], Theorem 1 associated to the measure p, 
which is the limiting spectral measure of our sequence (I3.24[) . in the sense of [S] . 
Thus one considers the spectral measure p\ associated to p.24p defined by 

rx{u -v)^ e[^ (^) ./^ {j)]=^ E Kx-(" " ^0 ■ ^^-^^^ 



Hence 

K 



k=l 

where dz stands for the Dirac measure at z e M^, |z| = 1. 

Since, by assumption (D) and the construction in Section 2, the measures p\ 
become equidistributed for X,K-^ 00, the limiting measure p is he normalized 
Lebesque measure on the unit circle and D is the constant associated to the RWM 
discussed in the Introduction; i.e. a — AnD. Recall p.l6p . p.27p and take say 
£2 = R~'^ ■ From the preceding, we obtain 

N^^+^Jx = (p + o(l))i?2 
NfJx = {v + o(l))i?2. (3.30) 

f2 



/T2 

and 
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In view of p.ip . p.2p . we obtain finally from the choice of £2, that 

Nf^{D + oil))X^. 

Recapitulating the preceding, the order in which the various parameters are chosen 
is 

R, £2, /3, £3, K, M, £4, £, £i,B{K, e). 
We proved the following 

Proposition 1. Assume E taken in a sequence such that (D) holds for E ^ 00 
and also, for some fixed 7 < condition 1(7, B(E)^ with B{E) ^-Z^ go. Let 

fE=J2 ^(^■^) 

or, more generally 

fs = a^e{x.^) with = a^, \a^ \ = 1. 

Then the number Ne of nodal domains of satisfies 

, Ne 

An > a 

E 

4. Arithmetic considerations 

We return to the assumptions (D) of equidistribution and (I) of independence. 
Recall that we assumed E of the form 

E = \{pl- (4.1) 

with Pa odd, Pa = l(mod4). 

Let TTa = |7r„|e*^" and write ^ = Ae^'^s ioi ^ e £ = £e, according to (HH]), pTTUl) . 

We start with a statistical discussion, considering a 'typical' integer E of the 
above form. 

A quantitative form of the required angular equidistribution is established in 
[EH] (see Theorem 1). 

Lemma 3. Given £ > 0, for almost all integers E considered above, one has a 
discrepancy bound 

AiE) < \£E\i\ogE)-''+' where K^^hg^. (4.2) 

Here A(_B) is defined as 

f3-a 



max 

O<a<0<27r 



2^ \£e\ - [#U e e [a,/3](mod27r)] . (4.3) 

The proof of this result depends on Kubilius' evaluation of the number of Gaussian 
primes in a sector and bounds on multiplicative functions. 
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Let us also recall that, on average, an integer E that is sum of 2 squares has 
X 5 log log -E prime factors, implying that \£e\ ~ y^logE. 

Next, we discuss the independence condition, again from a statistical perspective. 
The following statement follows from |B-Bj . Theorem 17 and Remark 15. 



Lemma 4. Given £ > 2, for most integers E of type (|4.ip . the number of non- 
degenerate relations (,i + ■ ■ ■ + = among elements of £e is at most 0{\£e\) 
for E ^ oo. More precisely, given any function Lp, most E , the 

number of non- degenerate solutions is bounded by if(\£E\)- 

Obviously, this implies property 1(7, B) for any given 7 > ^, for typical E taken 
large enough. 

Note that [B-Bj . Theorem 17 follows from the following statement, which in some 
sense is stronger. 
Denote 

^x.K ^{E^ \[p^ < X;pa = l(mod4),p„ > K}. (4.4) 
This set satisfies 

\nx,K\ - . (4.5) 

Vlog X vlog K 

Theorem 14 in |B-B| asserts then that for fixed i, 

lim lim — -\{E G fix K] £e admits a nondegenerate relation of length f\\ — 0. 

K^oD X^oo \\ix k\ 

(4.6) 

Lemmas [3] and 0] are clearly addressing the assumptions from our theorem in the 
statistical sense. Thus we can state 

Theorem 2. The conclusion from Proposition^^ holds for almost all E ^ 00 of 
the form pTi)) . 

Our next goal is a deterministic implementation. We start with the independence 
assumption. In certain cases, the desired information is provided by the deep work 
of Evertse-Schlickewei-Schmidt on additive relations in multiplicative subgroups of 
C* of bounded rank [E-S-Sj , which in turn depends on the subspace theorem. The 
result of ^E-S-S] states that any unit equation 

fliffi H h aege = 1 (4.7) 

with gi, . . . ,ge taken in a multiplicative group G over C of Q-rank r, has at most 
exp (c(£)(r + 1)) non-degenerate solutions. Here c{i) may be taken 

c{e) = {A£f^. (4.8) 

An immediate consequence is the following 
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Lemma 5. Let E = ria^iPa" '^^ above. Then the number of non-degenerate 
relations (|2.12p among elements from £e is bounded by 

exp(c(^-l)(2r+l))|£:£|. (4.9) 

While estimate (14.91) does not sufBce in general to conclude a condition 7(7, B), 
it does suffice provided r = o(log \£e\), i-e., recalling (|1.7p 

1 

- ^ logec cx). (4.10) 

Q — 1 

This is in particular the case if we fix the prime factors pi , . . . , of E and take 
their exponents Cq large enough. 

Remark 1 . It has been suggested that the true upper bound for the number of non- 
degenerate solutions of (|4.7I) may be subexponential in the rank, possibly bounded 

by 

exp {c{t)r^'^^^) for some < 1. (4.11) 
If this were true, 7(7, B) would hold with any 7 > -j, for all sufficiently large E. 

Remark 2. If we fixpi, . . . = l(mod4) and let E = ■ ■ ■ p'^^ with maxca — 00, 
condition (I) is certainly satisfied. Condition (D) amounts by (jl.lOp to equidistri- 
bution of the angular set 

r 

{2^ja0a;O<j„<e„} (modi). (4.12) 

a = l 

Hence 9^ ^ 27rQ. This is the case, since otherwise cos2b9a = 1 for some b € Z+, 
implying that cos26'q. is an algebraic integer. But since tgOa = |^ G Q, cos29a G Q 
and therefore cos2^?c( E I', 0^ £ (contradiction). 

Hence 

Theorem 3. Assume given pi,...,pr = l(mod4). Then the conclusion from 
Proposition [7] holds when E ranges in the set 

{P? ■ ■ -Pr'.ei, ...,ere Z+}. 

Remark 3. We may also state the following property, which results from (|4.7p . ()4.8|) 
and an easy adaptation of the proof of Lemm£[2l 

Theorem 4. Fix r and let E range in a sequence of integers of the form 

r 

E^Hp'^T (Pa = l(mod4)). 

For each E, let 
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be arbitrarily chosen, subject to the assumption that the probability measures 

PE=Y1 l«d'^A-e (A^ - E) (4.14) 

071 the unit circle, converge weak* to the normalized Lebesque measure on for 
i? — > oo. 

Denoting Ne the number of nodal domains of E , we have that 

Ne 

— > V for E ^ oo. (4.15) 



5. Further comments 

An alternative approach would consist in considering 'jets' of eigenfunctions. 
Thus given an eigenfunction of T^, introduce at each point x ^T'^ the scaled 
function 

^(y) = ^.(y) = /(a; + |). (5.1) 

The function ip may be e-approximated on [|y| < R\ by truncation of its Taylor 
expansion at order B = B{R, e), leading to a jet 

Jx = {D"ipx\y=i)}\a\<B- (5-2) 

Consider as a random vector in x e T^. Under assumptions (D) and (I), one may 
then show that the distribution of {Jx)x^T^ is approximatively the same as for the 
Gaussian random function with circular spectral measure and derive from this that 
^ V. This approach has the advantage of at least conceptually generalizing to 
real analytic compact manifolds M. Following [S , Section 2, one considers a map 
(assuming dimM ~ 2) 

= exp^ o4 : M2 ^ M, $,(0) = x 

with exp^ : T^M — > X the exponential map and : — > Tx{M) a linear Eu- 
clidean isometry. The function Lp^ is then defined by 

^x{y) = f{<^x{\''y)). (5.3) 

But we preferred to follow the procedure adopted earlier because it is more explicit 
and, in any case, we do not have examples at this point, other than the flat torus, 
where the RWM may be implemented deterministically. 

This discussion may however be of interest in the (arithmetic) hyperbolic case. 
(See |G-R-S) for some remarkable new results on nodal domains in this setting). 

Basically, the required behavior of the {Jx)x<em may here in some sense be seen 
as a far generalization of the Gaussian distribution conjecture of the eigenfunctions. 
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